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There is growing evidence that when magnetic reconnection occurs in high Lundquist number plasmas such 
as in the Solar Corona or the Earth’s Magnetosphere it does so within a fragmented, rather than a smooth 
current layer. Within the extent of these fragmented current regions the associated magnetic flux transfer 
and energy release occurs simultaneously in many different places. This investigation focusses on how best 
to quantify the rate at which reconnection occurs in such layers. An analytical theory is developed which 
describes the manner in which new connections form within fragmented current layers in the absence of 
magnetic nulls. It is shown that the collective rate at which new connections form can be characterized by 
two measures; a total rate which measures the true rate at which new connections are formed and a net rate 
which measures the net change of connection associated with the largest value of the integral of E\\ through 
all of the non-ideal regions. Two simple analytical models are presented which demonstrate how each should 
be applied and what they quantify. 
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I. INTRODUCTION 

The process of magnetic reconnection underpins our 
understanding of many astrophysical phenomena. Ex¬ 
amples include solar flares, geomagnetic storms and saw 
tooth crashes in tokamaks 1 2 . Yet a complete under¬ 
standing of this enigmatic plasma process remains illu¬ 
sive, despite decades of research. 

Fundamentally, magnetic reconnection is the process 
whereby excess energy in a magnetic field is liberated 
by the reorganization of a magnetic field’s connectivity 
in the from of plasma heating, bulk fluid motions and 
particle acceleration. Classically, this is envisioned to 
occur in a single well defined region of high electric cur¬ 
rent, within which non-ideal effects dominate and the 
plasma becomes decoupled from the magnetic fielcP®. 
However, in recent years the importance of instabilities 
which fragment reconnection regions has been more fully 
appreciated. In particular, in two dimensions high aspect 
ratio current sheets have been shown to be highly unsta¬ 
ble to tearing with the resulting dynamics dominated by 
the formation and ejection of magnetic islandJ^, whilst 
3D simulations have emphasized the importance of flux 
rope formation, braiding and the possible development of 
turbulenc^®33. Observations of plasma blobs and bursty 
radio emissions in the extended magnetic field beneath 
erupting CME’s as well as bursty signatures of recon¬ 
nection in the Earth’s magnetotail appear to somewhat 
corroborate this picture 11 14 . 

An important diagnostic of any reconnection scenario 
is the rate at which the process occurs. In two dimen¬ 
sions reconnection occurs only at X-points, with the rate 
of reconnection given simply by the electric field at this 
position. If the current layer is fragmented then the only 
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topologically stable situation is one in which only a single 
X-point resides at the boundary between the global flux 
domains. The reconnection rate is then the electric field 
measured at this “dominant” X-point (e.g. Wyper and 
PontinES). 

In three dimensions (3D) the picture is more complex. 
When reconnection involves 3D nulls, separatrix surfaces 
divide up the magnetic field into regions of differing con¬ 
nectivity. The rate of reconnection can then be defined 
as the flux transfer across these surfaced^, or past sepa¬ 
rators which sit at the intersection of different separatrix 
surfaced. If the non-ideal regions spanning the separa¬ 
trix surfaces are fragmented then consider ing flu x trans¬ 
fer across segments of a separatrix surfacd^® or along 
multiple separators 20 ! if they exist allows the reconnection 
rate to be quantified. Unlike 2D, where X-points other 
than the dominant X-point do not directly contribute to 
the reconnection rate (although they may indirectly af¬ 
fect it), in 3D reconnection across a separatrix surface in 
multiple places or at multiple separators all contribute 
towards the total rate of flux transfer between the main 
topological domains. This leads to the surprising result 
that in 3D two measures of reconnection may be used 
when reconnection occurs in fragmented current layers. 
One that measures the total rate at which flux is recon¬ 
nected (taking account of recursively reconnected mag¬ 
netic flux) and a net measure of the combined effects of 
each of the fragmented non-ideal regions. The former is 
the true reconnection rate for any problem, but the lat¬ 
ter may be of interest when the large scale effects of a 
reconnection site are being considered. 

Furthermore, in 3D reconnection may also occur in the 
absence of magnetic null points. In this case the lack 
separatrix surfaces against which reconnection can be de¬ 
fined requires a more general approach to the problem. 
The theory of General Magnetic Reconnection (GMR) 
encompasses reconnection across separatrices 1 ^^ as well 
as describing reconnection in situations without them. 
The theory of GMR has shown that for a single isolated 
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FIG. 1. Schematic of the difference in the evolutions of a 
and ^ seen by plasma elements on either side of the non-ideal 
region, D r . 


non-ideal region the rate of reconnection is given by the 
maximum off E\\dl on all field lines threading the non¬ 
ideal region 212 2J. However, the question remains as to 
how to measure reconnection in fragmented current lay¬ 
ers without the presence of separatrix surfaces or in situ¬ 
ations where separatrix surfaces are difficult to identify. 
The aim of this work is to extend the framework of GMR 
to quantify the reconnection process in this case. 

The paper is structured as follows. In Sec. II, we 
review the theory of GMR and introduce the relevant 
mathematical tools. Section III contrasts the manner in 
which new connections are created for single and multi¬ 
ple reconnection regions. Sections IV and V recap the 
derivation of the reconnection rate for an isolated region 
and then derive expressions for the reconnection rate in 
fragmented current layers. In particular, we show that as 
with reconnection involving null points a total and a net 
rate may be defined. The interpretation of each is then 
discussed. Section VI demonstrates the developed the¬ 
ory for two simple kinematic examples. Finally, Sec. VII 
summarises the new results and presents our conclusions. 


II. OVERVIEW OF GENERAL MAGNETIC 
RECONNECTION 

GMR is most readily developed within the framework 
of Euler potential^. A pair of Euler potentials (a and 
/3 say) are scalar functions which locally describe regions 
of non-vanishing magnetic field through the relation 

B = Va x V/T (1) 

As long as field lines are simply connected and only enter 
and leave through the boundaries of the region of inter¬ 
est once, a and f3 are single valued and can be used to 
label individual field lines, a and B are also flux coor¬ 
dinates and are related to the magnetic flux through a 
given surface via 


<f> = J J dad/3. (2) 

Coupled with an arc length (s) satisfying (B ■ V)s = B, 
any position within the volume of interest can be ex¬ 
pressed in (a, /3, s ) space. Within this formulation the 


electric field can be expressed as 

dot dB 

E —aF v/i + i v “- v *' < 3) 

where the quasi-potential (so named as it contains a 
time varying component) is related to the electrostatic 
potential B via 

^ = ( 4 ) 

when the magnetic vector potential is assumed to take 
the form A = aV/3. See Hesse and Schindler 22 for a 
discussion of the dependance of GMR on the choice of 
gauge taken for A. 

For maximum applicability a general form of Ohm’s 
law is assumed where the contributing non-ideal terms 
are grouped together into a single vector R such that 

E + v x B = R, (5) 


where R is assumed to be localized within a small region 
inside the domain of interest. By expanding R in covari¬ 
ant form and inserting it into Eqn. (J5| along with Eqns. 
0, § and © eventually leads to an expression giving 
the relative difference between the evolutions of a and (3 
that are locally 
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tpi and tp 2 are the quasi-potential functions on either side 
of the non-ideal region. Eqns. © and 0 show that 
plasma elements initially on the same field line threading 
a localized non-ideal region, D r measure a different evo¬ 
lution of a and j3 and so are not connected by the same 
field line at a later time. A sketch of this idea is shown 
in Fig. [l] The power of Eqns. 0 and 0 is that by 
considering only the relative difference in the evolutions 
of plasma elements, ideal flow components are removed, 
leaving only the components resulting in changes of field 
line connectivity and thus reconnection. If there is no 
variation in Z in a/3 -space then the evolutions of plasma 
elements are the same on both sides of D r . In this case 
plasma elements which begin on the same field line (and 
so initially have the same value of a and /3) are subject 
to the same change in a and ft and so will be found on 
the same field line at a later time. Therefore, a neces sary 
and sufficient condition for reconnection is that! 2 3 -^ 


V a jZ(a,l3)^ 0 , ( 10 ) 

i.e. that there be gradients in Z from one field line to 
another. 
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FIG. 2. Schematic showing the quasi-potential E mapped into flux coordinate space, (a) A single non-ideal region has one 
maxima of S. (b) Multiple non-ideal regions with multiple maxima and minima in S. Arrows show the direction along which 
new connections form. 


III. THE NATURE OF THE CONNECTIVITY CHANGE 

To understand the nature of any resulting connectivity 
change for a given problem it is useful to map the problem 
from 3D real space into flux coordinate (a/3) space. We 
will work in this space repeatedly throughout the rest of 
the paper. 

When the reconnection region is assumed to be local¬ 
ized within a single isolated region of E\\ the contours 
of E in flux coordinate space form closed loops. The 
Hamiltonian nature of Eqns. (|6| and 0 dictates that 
new connections be formed tangential to the contours of 
E. Thus, when E has only a single extrema these new 
connections will form in a circular manner. Figure 0a) 
shows a sketch of this concept where the green arrows 
indicate the direction along which new connections form. 

However, when a single reconnection region has an in¬ 
homogeneous E || or multiple reconnection regions exist 
within the region of interest then the mapping of E in flux 
coordinate space contains multiple maxima and minima, 
Fig. [2|b). In general we restrict ourselves here to sce¬ 


narios where the multiple E\\ regions still only make up 
a small fraction of the volume under consideration. This 
means that E approaches zero outside of a flux tube en¬ 
circling the multiple reconnection sites. The new connec¬ 
tions which form now do so along multiple closed loops 
embedded within a larger scale set of loops, Fig. 0b) 
(right panel). 

The way that this connection change is achieved de¬ 
pends upon the global constraints of the system under 
consideration. In general the formation of new connec¬ 
tions along these loops is a weighted combination of two 
extremes: steady state and purely time dependent con¬ 
nection change®: It is instructive to consider each in 
turn. 

In steady state the electric field is potential and the 
magnetic field remains fixed in time. Considering again 
the case when E has a single extrema, let us then assume 
that E = 0 on one side of the non-ideal region. The only 
way that new connections can form in the manner shown 
in Fig. |2ja), whilst also maintaining SB /dt = 0 is by 
inducing a circular plasma flow of the form shown in Fig. 
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FIG. 3. The two extremes of connectivity change. Steady state reconnection for a single non-ideal region (a) and multiple 
non-ideal regions (b). Purely time dependent reconnection for a single non-ideal region (c) and multiple regions (d). Orange 
denotes non-ideal regions where ^2 — 'ipi is positive and blue where it is negative. Solid arrows show the direction in which new 
connection form, whilst dashed arrows show induced plasma flows. 


[3ja). Ideal flows may be superposed on both sides, how¬ 
ever the connection change of the magnetic field within 
this ideal transporting flow will remain the same. Hornig 
and Priest^ considered one such example of this sce¬ 
nario. 

Extending this concept to multiple reconnection re¬ 
gions, each individual non-ideal region will behave lo¬ 
cally like the single reconnection region shown in Fig. 
[3ja). The key difference is that now a subset of field lines 
thread through multiple reconnection regions. Thus, cir¬ 
cular plasma flows are induced on field lines leaving a re¬ 
connection region which then feed into other secondary 
regions further along the same field line. Each secondary 
region superposes a circular plasma flow on to the flow 
pattern associated with the field lines which thread into 
it. In some cases this will enhance the induced flow at 
the exit of the patchy reconnection volume. In others 
it will act to reduce it. Figure [3jb) shows a conceptual 
sketch of this idea. Thus, steady state patchy reconnec¬ 
tion within a localized volume gives rise to an induced 
localized rotating flow with multiple internal vortices on 


field lines threading out of the reconnection volume. As 
with the single reconnection region any background ideal 
plasma flow may be superposed on to this non-ideal flow. 

In the opposite extreme of purely time dependent re¬ 
connection the electric field is assumed to be zero on both 
sides of the non-ideal region. This is particularly rele¬ 
vant to the Solar Corona, e.g. Priest and Forbes^. In 
this case new connections can only be formed by a time 
dependent change in the magnetic field within D r . The 
circular nature of this connection change in situations 
with a single extrema in E implies that helical magnetic 
fields are formed in the process, Fig. [3jc). When the vol¬ 
ume under consideration contains several reconnection 
sites, each helical region of field may contain a subset 
of field lines which threads into other helical reconnec¬ 
tion regions. Fig. id) depicts this idea. This shows 
that patchy time dependent reconnection can generate 
(or relax) braided magnetic fields which are t houg ht to 
be important in the context of coronal heating 27 28 . 

In any given 3D reconnection scenario a combination 
of both manners of connection change are likely to occur. 



5 



FIG. 4. The evolution of a flux surface 7 traversing the non-ideal region and passing through the maxima of E (viewed in 
flux coordinate space), (a) 7 (dashed line) over plotting contours of E for a single non-ideal region at t = t\. The green point 
denotes the maxima of S whilst the arrows show the direction of connection change, (b) the evolution of the 7 surface either 
side of a single non-ideal region at some later time, t = ti + At. (c) the associated swept out areas of flux (ST(t) and £#(£)). 


IV. QUANTIFYING RECONNECTION FOR A SINGLE 
RECONNECTION REGION 


If a given magnetic field contains separatrix surfaces 
and separator lines then these topological structures can 
be used as a reference against which the rate of flux trans¬ 
fer may be measured. For instance, when reconnection 
occurs along a separator the rate of reconnection is sim¬ 
ply given by the integral of E\\ along the separator line 17 . 
However, in the absence of such structures a more general 
theory is required. Hesse, Forbes, and Birn 25 developed 
such a theory for an isolated single reconnection region, 
D r , extending those of previous works 22 23 1 Using a simi¬ 
lar approach we now reproduce their results before gener¬ 
alizing the theory to quantify reconnection with multiple 
reconnection sites. 

Without an obvious reference surface against which to 
measure Hesse, Forbes, and Birn 25 considered an arbi¬ 
trary flux surface (i.e. a surface comprised of magnetic 
field lines) which intersects the single region of parallel 
electric field and contains the field line along which the 
integral of E\\ is maximal. When mapped into a /3-space 
this surface appears as a line, which they call the 7 line. 
Figure B a ) shows a sketch of this concept, where the 
contours depict the quasi-potential E. 

Now, this flux surface is comprised of field lines em¬ 
bedded in the ideal regions on either side of D r . Gener¬ 
ally, in each ideal region the evolution is comprised of a 
background ideal transporting component (which by def¬ 
inition is the same on both sides) and a non-ideal recon¬ 
necting component. Without loss of generality we now 
focus on the non-ideal component by fixing the evolution 
of field lines threading into the non-ideal region to zero, 
i.e. €1 = (da/dt\i,d/3/dt\i) = (0,0). This is equivalent to 
using a coordinate system which moves with the plasma 
on field lines entering the non-ideal region, allowing the 
connection change to be entirely characterized by the evo¬ 
lution of the field lines threading out of the non-ideal 


region which evolve according to €2 = (da/dt I 2 , d/3/dt I 2 ). 

If 7 is then defined at some arbitrary time (t = U), 
then at some later time (t == t\ + At) the differing evo¬ 
lution on either side of the non-ideal region splits 7 into 
two new flux surfaces. In af3 space these appear as two 
lines, shown in solid blue and dashed black in Fig. ( 4 ^b). 
Note that as e\ — (0,0) one of these lines is coincident 
with the original 7 line. The two new surfaces overlap 
at the edge of the non-ideal region and at £ maa , (where 
V^E = 0 ) since at these places e\ = €2 = ( 0 , 0 ). 

The magnetic flux reconnected up to this time is simply 
given by the flux bounded within one of the two flux tubes 
formed by these two new flux surfaces, denoted 5z,(£) and 
S R (t) (Fig. § c)). Each flux tube must have the same 
cross sectional area 


S(t) = S L (t) = S R {t), (11) 


due to the rotational nature of the connection change. In 
flux coordinate space this area is equal to the magnetic 
flux within each flux tube, recall the nature of the Euler 
potentials (Eqn. ©)• The rate of reconnection is then 
defined to be the rate at which S(t) (representing either 
Sl(£) or *Sr(£)) grows at t = G, 


^ = lim {- 7 - / dad/3\ = lim { / e ndsl , ( 12 ) 

dt t^XdtJs 7 t^ tl \f ds J y J 

where e = e\ on one side of S and e = €2 on the other, 
n is the outward normal of the boundary dS. As t —)> t\ 
the boundary of S collapses to become the section of the 
7 line on one side of the peak in E, referred to hereafter 
as 71 . The integral around the boundary of 5 at t = t\ 
then becomes the superposition of integrals along 71 , i.e. 
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t—yt 1 
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where l is the arc length along 71 . In coordinate space 
the local normal to 71 is given by 


“O-points”) of the divergence free field defining the di¬ 
rection of new connection formation: 


dp 3aA 

“ ~dl' ~dl) 


(14) 


dE dE 
62 = 1 ~dp’d^ 


(18) 


whilst Eqns. 0 and 0 give that 


e 2 


e 2, 


dE dZ\ 
d/3' da) ' 


(15) 


Combining Eqns. ( 12 ), (13), (14) and (15) then gives the 


reconnection rate as 


f (dE dp dE da \ 

dt J 7l \dp dl + da dl ) 



Thus, for an isolated region of E\\ with a single maximum 
of E the reconnection rate is given by the value of this 
maximum. This can be interpreted as the rate at which 
flux is transferred in one direction across any arbitrarily 
defined flux surface 7 which intersects the non-ideal re¬ 
gion and includes the field line upon which the maximum 
of E occurs. 


V. GENERALIZATION TO MULTIPLE 
RECONNECTION SITES 

When there are multiple reconnection sites or inhomo¬ 
geneity of E \| within a single site there is likely to be 
multiple peaks in E. We now aim to develop expressions 
which quantify the rate of reconnection in this case and 
explain their interpretations. 


A. Expressions for Reconnection rate 

As discussed in Sec. m when there are multiple peaks 
in E new connections are formed along multiple embed¬ 
ded closed paths in the aP plane. Near positive ex¬ 
trema (peaks) of E the direction that this new connec¬ 
tion formation takes is clockwise, whereas for negative 
extrema (troughs) it is anti-clockwise, Fig. 0b). The 
places where there is no connection change occur where 
V a ^S(a, P) = 0. These correspond to the field lines 
not threading into any non-ideal region (the neighboring 
ideal field) and special field lines along which the net dif¬ 
ference in connection change along their length is zero, 
i.e. field lines along which the induced connection change 
from multiple reconnection sites cancels out. These spe¬ 
cial field lines sit at the critical points (“X-points” and 


In terms of the quasi-potential the “O-points” correspond 
to peaks and troughs of E, whereas the “X-points” oc¬ 
cur at saddle points. Figure [5ja) shows a sketch of this 
idea, where the green and pink circles show the position 
of the “O-points” and “X-points” respectively. The key 
idea here is that just like X-points divide up two dimen¬ 
sional magnetic fields into distinct topological regions, 
so also the rotational formation of new connections de¬ 
scribed by €2 is partitioned into localized rotational re¬ 
gions (with “O-points” at their centers) by a series of “X- 
points”. The different topological regions of the €2 field 
are shown in different colors in Fig. [5] to better illustrate 
them. 

To construct expressions for the reconnection rate we 
begin in the same manner as Sec. [TV|and consider a flux 
surface bounded by the field line along which E max oc¬ 
curs and some other field line in the nearby ideal region. 
Figure [5|b) shows this surface as a dotted line in flux co¬ 
ordinate space. Now, the only topological regions of €2 
that are straddled by this surface are the region within 
which E max is situated (light blue) and the regions within 
which the surrounding outer loops of connect change oc¬ 
cur (white and orange regions). The connection change 
within these regions gives the net rate at which new con¬ 
nections form in a rotational manner around the field line 
on which occurs. The same analysis as Sec. [TV|may 
then be applied to this flux surface, quantifying this net 
rate as 


d$\ _ z 

dt )net~ 


-(/ ■ 

\Ja,P 


Ends 


(19) 


Thus, the maximum of the integral of E\\ across a frag¬ 
mented reconnection region measures a net rate of ro¬ 
tational connection change and neglects the connection 
change associated with other extrema of E. 

To quantify the true rate of reconnection requires that 
the connection change associated with each of these other 
extrema also be taken account of. This can be achieved 
by considering a flux surface for each additional extrema 
bounded on one side by the field line at which the ex¬ 
trema occurs and on the other by the field line situated 
at the nearest saddle point (s.p.) of E (corresponding 
to each additional “O-point” and its nearest “X-point” 
of the e 2 field). Figure |5|b) illustrates this idea with a 
series of dashed lines in flux coordinate space. For each 
of these additional flux surfaces one can also apply the 


same analysis as Sec. IV to give the rate of connection 
change across the surface as 


(-) = f 

\dt J local J 


= I c P d ^ 


dl 


— ^local extrema ^nearest s.p.* 


( 20 ) 
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FIG. 5. The evolution of bounded flux surfaces in a field with multiple reconnection regions (viewed in flux coordinate space), 
(a) Differently colored regions indicate the different topological regions associated with the connection change. Green circles 
denote maxima/minima and pink circles saddle points of S. (b) dotted line: a flux surface bounded by the the field line along 
which 5 is maximum and another in the nearby ideal region, used to calculate (d$>/dt) n et (see text). Dashed lines: flux surfaces 
bounded by field lines that have local extrema in S and field lines at the nearest associated saddle point, used to calculate 
(d&/dt)tot (see text), (c) the associated swept out areas of flux for each bounded flux surface. 




FIG. 6 . Examples of flux surfaces ( 7 ) chosen so that the extrema and saddle points of 5 are connected in a certain way (see 
text for details). Green and pink dots denote extrema and saddle points respectively. Arrows depict the direction of connection 
change. Dashed lines show sections of 7 forming chains of maxima and minima. Dotted lines show sections of 7 connecting 
these chains with each other or the background ideal field, (a)-(b) the same E profile with two different choices for 7 . (c) a 
choice for 7 in a more complex field, (d) a choice for 7 in a degenerate 5 field with a circular maxima (pink dashed line). 


The total reconnection rate associated with all of the 
non-ideal regions is then given by the sum of the local 
connection change occurring around each additional ex¬ 
trema in addition to the net rotational connection change 
occurring around the largest extrema, i.e. 

) = |^max| T ^ ^ l^local extremal ^adjacent s.p.,i| •> 
tot 1 

( 21 ) 

where absolute values have been used to account for when 
E has both maxima and minima. 

The above shows that when reconnection occurs in 
fragmented reconnection regions the rate of reconnection 
can be quantified by two different measurements. The 
first, {d<&/dt)tot describes the true rate at which new con¬ 
nections are formed collectively by the fragmented layer. 
The second, (d&/dt) net gives the net rate of flux trans¬ 
fer associated with the maximal value of E on field lines 



crossing the volume. Both measurements are equivalent 
when there is only one peak in E. 


B. Interpretation in Terms of Flux Transferred Across a 
Single Large Flux Surface 

The above analysis shows that reconnection in frag¬ 
mented current layers can be considered as representing 
the rate at which magnetic flux is reconnected across mul¬ 
tiple bounded flux surfaces. We now show that provided 
E is smooth and continuous, (d<$>/dt) to t and (d$/dt) ne t 
can also be interpreted as the total and net rate of flux 
transfer across a large scale flux surface ( 7 ) spanning the 
entire reconnection volume. 

From the above analysis one would expect that such 
a flux surface must contain the field lines along which 
each of the extrema and saddle points of the E profile 
occur. However, the order in which each extrema and 





























FIG. 7. Interpretation of reconnection rate in terms of a global flux surface, (a) Dashed line shows a flux surface ( 7 ) which 
traverses all of the non-ideal regions and includes the field lines at which the extrema and saddle points of S are found, (b) 
variation of 5 as a function of arc length l along the 7 line in flux coordinate space. Li and Ri denote the direction of flux 
transfer across 7 (where i G [1,4]) and the X’s and O’s show the positions of the saddle points and extrema of 5 respectively, 
(c) the evolution of the 7 surface either side of the multiple non-ideal regions and (d) the associated swept out areas of flux, 
Sl,z( t) and SR,i(t) corresponding to the different sections of flux transfer across 7 (Li and Ri). 


saddle points are connected by 7 is crucial. In partic¬ 
ular, 7 must be defined such that extrema of the same 
type (maxima or minima) are connected via any adjoin¬ 
ing saddle points forming chains. The end of a chain of 
maxima can be connected with the end of a chain of min¬ 
ima if the connection is from the maxima to minima, or 
saddle point to saddle point. Alternatively, the ends of 
chains of maxima or minima may instead be connected 
with the surrounding ideal field. Figure J 6 ^a)-(c) shows 
three examples. In the degenerate case of when a local 
maxima of 5 forms a ring, any two points on the ring 
may be chosen in place of two saddle points, Fig. §<*)• 

We point out that the selection of this flux surface 
is not unique and differs depending upon how different 
chains of maxima and minima are connected. Figure 
J 6 ^a)-(b) illustrates this idea by the differing dotted sec¬ 
tions of 7 . 

The dashed line depicted in Fig. [7|^a) shows a flux 
surface which connects the extrema and saddle points of 


the previously considered £ profile in the way described 
above way at some arbitrary time (t = t\). By choosing 
7 in this way the flux transfer between the critical points 
of E alternates along 7 . This is shown in Fig. (7^b) which 
depicts the variation of E along 7 as a function of the arc 
length (i). The positions of the extrema are shown with 
“O’s” and the saddle points with “X’s”, whilst the direc¬ 
tion of flux transfer across 7 between them is indicated 
by or “Li”, where i G [1,4]. 

At some later time (t = t% + At) the differing field 
evolutions on either side of the multiple non-ideal regions 
forms a chain of flux tubes with cross sectional areas of 
SL,i(t) or SRj(t) associated with Li and Ri respectively. 
Note, that the rotational nature of e 2 means that the 
sum of each set of area elements must be the same, i.e. 
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FIG. 8 . Reconnection of flux across a specific flux surface, (a) 
dashed line depicts 7 in flux coordinate space. The contours 
depict E and the colors denote the topological regions associ¬ 
ated with 62 . (b) E as a function of arc length l along 7 . Li 
and Ri (where i £ [1,3]) denote the direction of flux transfer 
across 7 . 


Now, if we compare the area segments swept out by the 
series of bounded flux regions discussed earlier (Fig. J5^c)) 
to those generated by this continuous surface (Fig. |7fb)) 
we find that they match SL,i{t ). This shows that the 
total rate of reconnection can be interpreted as the rate 
of growth as t —)> t\ of the collective area S(t) associated 
with flux swept in the same direction (all to the left or 
all to the right) across 7 , i.e. 



A similar conclusion can also be drawn for other choices 
of 7 connecting the chains of maxima and minima. 

Similarly, the net rate can be interpreted as the rate of 
growth as t t\ of the difference in the areas associated 
with flux swept in one (or the other) direction on one 
side of l^ maxi be. 



net 



(24) 


where 


S d {t) = 


?(/... 


dad/3 > — ^ 


' SR,k 


dad/3 > , (25) 


and j and k sum over the area segments formed along 
the portion of 7 on one (or other) side of E maa ,. 

It should be noted that the existence of such a large 
scale 7 surface is not necessary for the application of 
Eqns. (19) and (21), and indeed if E is sufficiently com¬ 


plex or contains discontinuities such a surface may not be 
definable. However, we have shown that at least when E 
is smooth and relatively simple the intuitive idea that the 
reconnection rate should measure the rate at which flux 
is reconnected across some large scale flux surface (akin 
to that of a true separatrix when reconnection occurs 
between distinct topological regions) still holds. 


C. Quantifying Reconnection Across and Arbitrary Flux 
Surface 


Finally, we now consider the case where rather than 
wanting to know the true rate of reconnection, one is in¬ 
terested in knowing the rate at which flux is reconnected 
past a particular flux surface. An example of such a sur¬ 
face would be one associated with an observed flare rib¬ 
bon on the photosphere. Another would be if the global 
topology is such that field lines from a separatrix sur¬ 
face or surfaces pass through the domain of interest and 
one wishes to know the rate of flux transfer between two 
different topological domains. 

Equations ([l9|) and (21) are easily generalized to this 


scenario. Consider some arbitrary flux surface spanning 
a fragmented reconnection region with multiple peaks in 
5, Fig. § a). Along the length of 7 a number of local 
maxima and minima of E occur. Between each of these 
local extrema flux is transferred in one or other direction 
depending upon the sign of the gradient of S(/), Fig. 
| 8 jb). In analogy to the previous sections the net rate at 
which flux is transferred across this surface is given by 



where the subscript 7 denotes measurement of each quan¬ 
tity along the 7 line. Similarly the total rate of flux 
transfer across this particular flux surface is given by 


— |^max,7 |T ^ ^ l^local max,i ^adjacent min,z|-y 
tot, 7 1 

(27 

Depending upon the path take by the 7 line as it crosses 
E in flux coordinate space the value of (dE/dt)tot^c an b> e 
greater or less than the value measured by Eqn. pi) ). For 
instance if 7 is chosen so that it crosses E many times, 
then it would be likely that (dE/dt)tot , 7 > (dE/ dt) tot . 
However, by definition the net rate of transfer will at 
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FIG. 9 . Iso-surfaces at 10 % of the maximum of |R|, showing 
the three localized non-ideal regions at t = 0 in both models. 
In red are a selection of field lines plotted from footpoints 
along (x, z) = (0, 2). 


most be the same as the net rate of rotational connection 
change around the field line with E = E rnax , so that 

(6^/^)net,7 {d^/dt) ne i. 


TABLE I. Model Parameters 


i 

ji 

lx,i 

ly,i 


%0, i 

2 /o,i 

zo,i 

1 

-0.1 

0.2 

0.2 

0.2 

O 

o 

0.0 

0.0 

2 

-0.1 

0.1 

0.1 

0.1 

0.0 

0.35 

-1.0 

3 

-0.1 

0.1 

0.1 

0.1 

o 

o 

- 0.35 

-1.0 


VI. EXAMPLE: MULTIPLE RECONNECTION SITES IN 
A STRAIGHT FIELD 

To illustrate the theory we now present two simple 
kinematic models of an idealized fragmented current 
layer. 

Starting with an initial magnetic field (at t = 0) of the 
form 


B = 5 0 z, (28) 

we assume some non-ideal process occurs to produce mul¬ 
tiple non-ideal regions such that 

n 

i =0 

(29) 

where (l X}i , l V}i , l zA ), (x 0 ,i, yo,i, z o,i) and i* control the di¬ 
mensions, position and the strength respectively of each 
non-ideal region. We choose three non-ideal regions 
(n m 3), one larger central region and two smaller iden¬ 
tical offset regions, see Fig. [9j The chosen parameter 




y 


FIG. 10 . (a) contours of E(x, y) mapped on to flux coordinate 
space, (b) E along the line x = 0 , passing through the five 
critical points. 


values are given in Table [T| Depending upon the con¬ 
straints placed upon the system, reconnection solutions 
describing purely time dependent, steady state or a com¬ 
bination of both scenarios can be constructed. In what 
follows we will consider the two extreme cases and verify 
in each case the validity of the Eqns. (19) and ( 21 ). 


A. Time Dependent Reconnection 

In this extreme we impose that the sections of field 
lines threading into and out of the non-ideal region are 
held fixed such that the electric field vanishes on each side 
of the non-ideal region. This is equivalent to assuming 
that the plasma velocity v = 0 everywhere. Ohm’s law 
then gives directly that E = R, i.e. 

3 

E = ^ j. e -( x ~ x o,i) 2 / l l,i-(y-yo,i) 2 / l 2 y :i -( z - z o,i) 2 / l l,i 2. 

i =0 

(30) 

Faraday’s law, d~B/dt = V x E then dictates that at later 
times the magnetic field evolves such that 

3 

B = Bq Z T ^ ^ V X A i, fluxring •> (31) 

i =0 

where i sums over each non-ideal region and 

A , 7 . -fj. p -( X - X 0,i) 2 / l l,i-(y-y0,i) 2 / l l,i-( Z - Z 0,i) 2 / l2 z,iy 

■^i^jluxring — Zi - 

(32) 

At t = 0 the magnetic field is initially straight, but as 
time progresses each flux ring introduces an ever increas¬ 
ing twist to the field. Note that in this simple example 
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we are only considering small periods in time, t. At t = 0 
the straight magnetic field can be described with the two 
Euler potentials a = x and (3 = y. Since each non-ideal 
region is negligibly strong in the vicinity of the others 
E(a,/3) = E(x,y) can be constructed from the super¬ 
positions of E\\dl = E z dz across each region 
giving 


ri 

F(x, y) = JtL e -( x - x o ,i) /i x ,i~(y-yo,i) / l v,i-\-f(x,y) 


i =o 


(33) 

where f(x,y) is an arbitrary function. In what follows 
we will trace field lines from z = +2 to z = —2 so for 
convenience we set E(x, y) = 0 at 2 = +00 to give 

E(x,y) = ^-e-^- X0 - i)2/l2 ^ v - V0 ^ 2/l v^ (34) 


n 

1=0 ’ 

Figure [lO^ shows a contour plot of S at t = 0 mapped 
on to the xy- plane. The profile contains three distinct 
peaks (O-points) with two saddle points (X-points) be¬ 
tween them. By symmetry the X-points and O-points of 
€2 — 61 = €2 lie along x = 0 , so we choose this as our 7 
line. 

The variation of the quasi-potential along this line 
E(x = 0,?/) is shown in Fig. 10 3 . The peaks occur at 
y = a, c and e, with the saddle points located at y = b 
and d. Applying Eqn. (21) gives the total reconnection 
rate of this system as 



FIG. 11 . (a) connectivity map at t — 1 for the time dependent 
model. Black show field lines with starting points below x = 0 
and white those with starting points above, (b) connectivity 
plot of the field lines to have changed connection between 
t — 0 and t — 1 . Black regions have moved from x < 0 to 
x > 0 , white have moved from x > 0 to x < 0 and grey regions 
have stayed the same. 


where N is the number of field lines under consideration. 
The rate of reconnection is then estimated as 


d§ _ A<1> 
dt At 


(38) 


) = Ew + £ |S local extrema;i - ~ adjacent s ,. ;i | 

/ tot i 

— 'c T |^a ^b\ T |^e ^d\ 

« 0.04711, (35) 

with a net rate of flux transfer given by 

= Zmax ~ 0.03545. (36) 

net 

In this extreme, these values represent the total and net 
rate respectively at which magnetic field is generated nor¬ 
mal to the 7 surface collectively by the non-ideal regions. 

We now go on to verify these values by comparing them 
with values obtained numerically from a flux counting 
procedure, explained below. A large number of field lines 
were traced from a grid on 2: = 2 as far as z = —2. At 
both positions the magnetic field has reached its asymp¬ 
totic value of B = B 0 z. This is done for the field at 
some time, t = t\ and some later time t = t\ + At. The 
amount of flux transfer (A<F) in this period is obtained by 
comparing the final positions (on z = — 2 ) at both times 
and summing the number of field lines to have crossed 
the 7 line, weighted by their area element on the starting 
grid and the field strength perpendicular to the surface 
of starting points, i.e. 

A<F = ^ Ax Ay, (37) 

N 




To obtain (d$/dt) to t all field lines found to have crossed 7 
in At are counted and the value halved so as not to double 
count the flux transfer (recall that the connection change 
is circular and so will cross the 7 line twice). (d$/dt) net 
is approximated by counting only the net transfer across 
a half segment of the 7 line. 

The mapping of field lines on 2 ; = —2 at f = 1, color 
coded according to whether they start above or below 
7 on the other side of the non-ideal region (z = 2) is 
shown in Fig. [llja). Figure [Tljb) shows the regions 
within which field lines have changed connectivity com¬ 
pared with the mapping at t = 0. White areas depict 
where flux has reconnected across 7 from x < 0 to x > 0, 
and black regions where flux has been reconnected in the 
other direction. Grey shows regions where field lines have 
not crossed 7 . Figure [T2| shows a 3D visualization of the 
field at t = 1, were the iso-contours depict the shape 
and position of each non-ideal region. Applying the flux 
counting procedure we obtain that 


d$ 

dt 


' d® 

°- 04857 ’ l* 


: 0.03607 


(39) 


for a grid of 400 2 starting points. Aside from a small 
variation due to the discrete nature of the method, these 
results agree closely with the value obtained by applying 
Eqns. (friT) and (HI)). 


Lastly, consider now the instantaneous reconnection 
rate at the later time (t = 1). At t = 1 each non-ideal 
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FIG. 12. Iso-surfaces at 10 % of the maximum of |R|, show¬ 
ing the three localized non-ideal regions at t = 1 in the time 
dependent model. In red are a selection of field lines plotted 
from footpoints along (x,z) = (0,2), demonstrating the in¬ 
jection of twist into the field and the overlap of the field line 
mappings. 



X 


0.4 



(b) 


- 1.0 - 0.5 0.0 0.5 1.0 

y 


FIG. 14 . (a) E(x, 2/) in the steady state example, (b) The 
induced perpendicular plasma flow at z = — 2. 



giving that 


= ~ J R-B/|S|ds + 0 o 
= — J E\\ds + 4 >o 


= E + 0o 


(41) 


For illustration we set 0o = 0 which removes background 
ideal motions. Thus, 


FIG. 13 . The quasi-potential calculated numerically at t — 1 . 


E = -VS. 


(42) 


region now adds a non-zero twist to the field line map¬ 
ping. The overlapping nature of the mappings distorts 
the shape of S and therefore the positions of the extrema 
and saddle points, Fig. [13] As a result the conceptual 
flux surface 7 against which reconnection rate is being 
measured by Eqn. (21) moves to pass through these 
points at this later time. 


This electric field differs from R, with a non-zero part 
outside of the non-ideal region which induces a perpen¬ 
dicular plasma flow of the form 


V _L 


E x B 

B 2 


(43) 


The magnetic field in this case remains straight for all 
time, and the quasi-potential is simply the same as the 
time dependent case at t = 0 , i.e. 


B. Steady Sate Reconnection 

For comparison we now consider the opposite extreme 
of steady state reconnection for the same initial magnetic 
field and non-ideal term (R). In steady state the electric 
field can be expressed in the form of a potential 

(40) 


E{x,y) = V^J2 (44) 

i =0 

Figure pA] shows the induced plasma flows on one side 
of the reconnection regions when the electric field is as¬ 
sumed to be zero at z = 2. The generated flux trans¬ 
porting flows follow the contours of the quasi-potential, 
producing three overlapping vortices. As the contours of 
S now form the stream lines of the perpendicular plasma 
flow, the zeros in the flow pattern are co-located with 


E = —V0 = —v xB + R, 
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FIG. 15 . (a) velocity perpendicular to the 7 line (x — 0 ) in 
the steady state example, (b) variation of S along the 7 line. 
Note that the zeros in the velocity field correspond to peaks 
or troughs of S. 


the peaks and saddle points in 5, Fig. [l5j As the quasi¬ 
potential £ is the same as the time dependent scenario 
at t = 0 the two measures of reconnection rate are then 
also 


« 0.04711 & (^) « 0.03545, (45) 

tot \dt J net 

where 7 can be chosen to lie along x = 0. In this extreme 
these quantities are measures of the total and net rate at 
which flux is swept past x = 0 by the induced plasma 
flow on one side of the collective non-ideal regions, i.e. 




/ oo 

B 0 \vj_(x = 0,z = -2)| +/-dy, 

-OO 

/»00 

= / B 0 v ± (x = 0,z = -2)dy, 

J 0 


(46) 

(47) 


where |..| + /_ denotes integration over either the positive 
or negative values only. An approximate expression for 
this flux transporting flow evaluated on 7 (x = 0 ) at 
z = —2 is 


v_L, 7 = v x (x = 0,z = -2) = EyB z /B 2 « ~^B 0 /Bq 

V (48) 

which when substituted into Eqn. ( [46] ) and integrated 
over the regions of negative velocity leads to 

(jP) tot = (Se “ “ d) + (Sc “ Eb) + (Sa “ 0) ’ 

= S c + |E a — Eb| + |S e — Sd|. (49) 


Note that integrating over the positive value gives the 
same result. Substituting the above expression for v_l ?7 
into Eqn. (47) then also gives that 


d$\ _ „ 

^)net~^ C 


(50) 


Eqns. ( [49] ) and ( [50] ) are simply Eqns. ( [ 21 ] ) and ( [19] ) 
applied to this particular E profile. 

Thus, we have verified the two rates of reconnection 
for our idealized fragmented reconnection region in each 
of the two extreme cases of steady state and purely time 
dependent reconnection and by extension the continuum 
of cases in-between. 


VII. DISCUSSION AND CONCLUSIONS 

The aim of this paper was to extend the theory of 
General Magnetic Reconnection to situations with frag¬ 
mented current layers within a localized volume. We 
considered the manner in which new connections may 
be formed, derived expressions for the rate at which this 
occurs and verified these expressions with two simple ex¬ 
amples. 

In terms of facilitating the formation of new connec¬ 
tions we showed that in the extreme of steady state re¬ 
connection a large scale rotational non-ideal flow with 
internal vortices is produced, whilst purely time depen¬ 
dent reconnection leads to spontaneously braided mag¬ 
netic fields. However, it should be emphasized that the 
reverse is also true. That is, the existence of non-ideal re¬ 
gions is guaranteed by the right evolution of the magnetic 
field (given the necessary non-ideal plasma conditions). 
In particular, if a magnetic field is initially braided with 
the field lines entering and leaving the volume held fixed, 
then multiple current layers must form to remove this 
braiding. This second scenario is readily observed by nu¬ 
merical experiments examining the non-ideal relaxation 
of braided magnetic fields (e.g. Pontin et alP Rappazzo 
and Parker^). 

By considering the closed paths along which these new 
connections formed we also showed that when current 
layers are fragmented two rates of reconnection can be de¬ 
fined which describe the process. (d&/dt) to t which mea¬ 
sures the true rate at which new connections are formed 
collectively by the multiple non-ideal regions and a sec¬ 
ond, (d$/dt) net measuring the net rate at which changes 
in the global field occurs. When applied to a single re¬ 
connection region both rates are equal. 

We chose to define (d<$>/dt) to t such that it measures the 
total rate at which flux is locally and globally cycled when 
viewed in flux coordinate space. This requires evaluating 
the quasi-potential at the saddle points of E as well as 
the extrema. We chose this rather than a simple sum 
over each extrema as summing over only the extrema 
overestimates the rate flux is cycled (although if each 
non-ideal region has little overlap this may give a close 
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approximation, e.g. Pontin et al. 13 ). This occurs as 
each extrema taken on its own measures the net rate of 
transfer of flux between itself and the background ideal 
field. Therefore, summing over all extrema double counts 
the flux being cycled around outer loops, such as those 
depicted in orange and yellow in Fig. [5] By involving 
the quasi-potential measured at the saddle points, this 
double counting is avoided. 

It is also worth emphasizing that our total reconnec¬ 
tion rate (d^/dt) to t does not measure the sum of the 
reconnection rates of each individual reconnection region 
within the volume. The only way that this could be quan¬ 
tified would be to consider the local quasi-potential drop 
across each non-ideal region in turn. However, each re¬ 
gion would have to be surrounded by ideal magnetic field 
for this to be meaningful. In fragmented current layers 
this is rarely the case as different current sheets partially 
overlap when merging or breaking apart. Considering the 
collective behavior as we have done here is the only way 
to properly quantify such a system. 

Given that we have introduced two different rates to 
describe this collective behavior, which should be used 
to characterise a given reconnection process? It depends 
upon what is most of interest for the problem at hand. 
For instance, if one is considering the scaling of energy re¬ 
lease compared with reconnection rate then the total rate 
is the better choice. It would also be the more relevant 
choice in situations where the rate at which flux is swept 
up by a fragmented reconnection region is of interest, as 
is thought to be related to photospheric brightening in 
solar flares (e.g. Hesse, Forbes, and Birn^, Qiu et a/.®). 
However, the net rate may be more useful when the mul¬ 
tiple reconnection regions are fluctuating and transient 
(as occurs during an increasing turbulent evolution of 
the magnetic field) and there are some simple large scale 
symmetries against which flux transfer is wished to be 
know (e.g. Wendel et Kowal et al 32 ). 

Ultimately the non-ideal physics associated with the 
plasma, any gradients in the mapping of the magnetic 
field and the way in which excess magnetic energy is built 
up will dictate where non-ideal regions form and if they 
subsequently fragment. The present analysis serves as 
a way of interpreting how the subsequent reconnection 
proceeds and how best to quantify it. 
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